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MICROWAVE  MEAT  ROASTING  - A COMPUTER 
ANALYSIS  FOR  CYLINDRICAL  ROASTS 


I.  INTRODUCTION 

*n)e  research  method  and  results  of  a program  in  microwave  meat  roasting  during 
1974  and  1975  has  been  reported  in  the  Journaj  of  Microwave  Power.'  The  key  element 
of  this  research  effort  was  the  development  of  a computer  model  of  microwave  meat 
roasting.  'Hiis  computer  model  consisted  basically  of  the  numerical  solution  of  the  heat 
conduction  equation  applied  to  the  node  network  of  a cylinder;  the  output  of  this  program 
consists  of  node  temperatures  for  each  time  increment  throughout  the  cooking  period. 
The  solution  of  the  temperature  distribution  throughout  a solid  cylinder  has  been 
previously  obtained  in  closed  form  for  various  initial  conditions  and  heating  schemes.^ 
The  complication  added  by  microwave  heating,  however,  precludes  the  use  of  such  closed 
form  solutions.  Also  complicating  the  issue  is  the  evaporative  cooling  at  the  meat  surface 
which  is  both  time  and  temperature  dependent.  This  report  describes  the  details  and 
assumptions  associated  with  the  computational  procedure,  as  well  as  detailed  material 
compiled  from  references  but  not  listed  in  reference  1.  This  report  is  intended  as  reference 
for  users  of  the  ROAST  computer  program. 

II.  HEAT  TRANSFER  EQUATIONS 

It  is  assumed  that  microwave,  convective,  and  radiant  heating,  and  evaporative  cooling, 
are  uniformly  distributed  over  the  cylindrical  roast  surface.  It  is  also  assumed  the  meat 

is  homogeneous  and  all  properties  are 
isotropic.  Symmetry  about  the  r and  z axes 
is  therefore  assumed  and  temperatures  in  the 
thin  wedge  (shaded)  are  representative  of 
temperatures  throughout  the  entire  roast. 
Excluding  terms  involving  6,  the  general  heat 
conduction  equation  for  a cylinder  is: 

pc  dT 

dr^  r dz*  dz^  k k dt  (1) 


' Nykvist,  W.  E.  and  Decareau,  R.  V.,  "Microwave  Meat  Roasting,"  J.  Microwave  Power 
11  (1),  1976 

’Schneider,  P.  J.,  Conduction  Heat  Transfer.  Addison-Wesley,  Reading,  MA  1965 
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T « temperature, 
t * Time,  s 

q « rate  of  heat  addition,  cal/t  cm’  . 

k " thermal  conductivity,  cal/cm  s°C 

p » densitv,  g/cm’ 

c * specific  heat,  cal/g^C 


Since  die  dielectric  properties  of  meat  ere  temperature  dependent,  the  penatration 
of  microwave  radiation  varies  with  temperature,  greatty  complicating  en  analytical  solution 
of  equation  (1).  This  equation,  therefore,  is  approximated  by  a finita  difference  technique. 

A network  of  nodes  Is  imposed  over  the 
center  plane  of  die  thin  wedge  and  die 
equation  is  rewrittan  using  finite  difference 
approximations  for  the  differentials.  Using 
• {m,n+l)  ^ twodimensional  node  scheme  shown, 

^ equation  (1)  in  finite  difference  form 

becomes 


T 


^m+1,n  ^m*1,n  ~ 2^m,n  ^ ^m+1,n“'^m*1,n  '‘'^m.n+l  ^m,n*1  ~ ^^m,n 

r 


Ar’ 


2Ar 


+ q . pc  r(^rn,n^new  ” ^m,n"| 

k ITL  At  J 


12) 


Rearranging,  letting  q * p/V  and  solving  for  tha  new  node  temperature,  9'''®’* 


^^m,n^new 


'*'m+l,n  ^m*l,n  ^m+l,n  ^m*1,n  ^m,n-l 

pc,  Ar*  2rAr  Az’ 


+ P + T„„f^®-  2 

^ m,n  * 


kV 


kAt  Ar’ 


2 

Az’ 


(3) 


where  p ~ microwave  power  absorbed  by  node  m,n  (cal/s) 


V * volume  araociated  with  node  m/i  (cm’) 


Hotman^  provides  a detailed  description  of  finite  difference  approximations. 


A closer  view  of  the  thin  wedge  shows  each  node  to  have  a volume  associated  with 
it,  as  shown  in  Figure  1. 


(1.1)  (1.2)  f1.S) 


Figure  1.  Roast  Wedge  Showing  Typical  Node  and  Associated  Volume 

Figure  1 has  a 5-by-5  network  of  rKXles;  in  the  ROAST  computer  program  nodes 
are  usually  spaced  every  0.5  cm,  giving  typically  a 13^y-25  node  network.  Once  the 
roast  dimensions  and  the  desired  number  of  radial  nodes  is  known,  the  “end  areas"  denoted 
by  A(m)  are  calculated  by  assuming  the  arc  length  at  the  roast  surface  to  be  equal  to 
Ar.  Thus,  associated  with  each  node  are  arc  lengths  S(m)  and  Slm-I),  and  area  Aim). 


^Holman,  J.  P„  He$t  Transfer,  McGraw-Hill,  New  York,  1968 

j 
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There  are  two  types  of  boundary  conditions  for  the  wedge,  the  inside  or  noeat-to-noeat 
interface  and  the  outside  or  meat  surface-to-air  interface.  The  inside  boundary  conditions 
ere  handled  in  subroutine  INTER  as  special  cases  of  equation  (3).  There  are  three  special 
cases  for  insida  boundary  conditions; 

Casa  1 m » n B 1 


This  case  involves  the  wedge  innermost  comer  and  involves  heat  transfer  from  only 
two  adjacent  nodes.  To  take  into  account  symmetry  at  tha  two  wedge  edges,  "ghost" 

ghost  assumed  to  be  outside  the  wedge. 

These  nodes  are  assumed  to  be  at  the  same 
temperature  as  the  real  nodes  opposite.  The 
form  of  equation  (3)  then  becomes 


2T||in+!  + 


P 

kV 


'm,nJ 


Case  2 m - 1,  1 < n < MHL 

•(m+1,n)  ghost 


This  case  covers  nodes  along  the  top  edge  of 
the  wedge.  Using  only  ona  ghost  node, 
equation  (3)  for  this  special  case  becomes 


' m,n+l 


rm,n-1 


Az* 


P 

kV 


Vn/£1__2  _ jj 

IjcAt  Ar^  A2*J 


15) 


Cbm  3 n*l,  l<m<  MR 


(T 


wedge.  The  appropriate  form  of 


*^’m+1,n  ~ ^Tm  n+I  + P + 

kV 


2rAr 


A 


(6) 


Outside  boundary  conditions  involve  heat  transfer  at  the  meat  surface.  These 
conditions  are  handled  by  making  a heat  balance  about  each  surface  node  instead  of 
modifying  equation  (3).  The  three  basic  formulas  used  are 

At 


(a)  heat  conducted  = k (conduction  area] 

(path  length] 

(b)  heat  convected  * h (surface  area]  ITm,n-Tal 

(c)  increase  in  internal  energy  > pc  (node  volume) 


_ 


^^m,n)new  ^m,n 


At 


-J 


A special  subroutine,  BNDRY,  handles  the  following  five  special  cases; 

Case  4 m ■ 1,  n * MHL 


(•T\0-1| 


(nW'l 


This  case  covers  tha  node  at  the  upper  right 
corner  of  the  wedge.  No  ghost  node  need 
be  considered  in  this  case  as  a direct  heat 
belence  is  being  carried  out.  This  node  can 
receive  radiant  and  convective  energy  and  can 
low  heat  by  evaporative  cooling  at  tha  outside 
surface  designated  by  A(m).  The  heat 
balance  is 
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Kestfrom 
to  (m,n) 


heat  from 

(m+l,n)  to  (m,n)  + convected  heat  + 


kA(m) 


|_Az_  2_  Ar 


radiant 

microwave 

evaporative 

heat 

+ heat  - 

haat 

= increase  in  internal  energy 

drAlm) 

+ p - 

hgA(m) 

“ ^^f<Tm.n>new 

2 L At 

This  equation  is  solved  for  (T^  „>^cw 


^'*'m,n^r 


pcA(m)  A z 


Vn  l ^ S(m)Az  T^->-l,n  * T,  + q 

Az  2Ar  k 

* b 


qfAlm) 


+ p - h^(m)  + Tn,  „ JpcA(m)  Az  - A(m) 
k k \“  2A  tk  ~ "aT 


— S(m)  Az  — hA(m) 
2Ar  ~k 


(m,n-1)^ 


Cm  5 1 < m < MR,  t?  * MHL, 

•{m-l,n) 

It  4--» Aim) 


This  case  covers  the  nodes  on  the  right  side 
of  the  wedge.  Again,  tha  surface  with  area 
Aim)-  is  subject  to  radiant,  convective  and 
evaporative  heat  transfar. 


— Sim) 
•|m+l,n) 
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Writing  tha  heat  balance  and  solving  for  (T^,n^new  gives 


tTni,nV 


rA(in)  + S(»n)A2  + S(m'1)Az 

pcA(m)Az[_Az  2Ar  2Ar 

hA(m)Ta  + QrAlm)  + p - h^|m)  + T,n  n^pcA(m)Az 
k k k 1 2Atk 


pcA(m)Az 


- Aim)  - S(m-1)Az  - hA(m) 


Case  9 m MR,  m • 1 


•Un-l,nl 


This  case  covars  the  lower  left  comer  of  the 
wedge.  Hara  tha  curved  surface  with  area 
S(m)Az/2  is  the  outer  heat  exchange  surface. 
The  heat  balance  becomes 


•(m,n+1) 


*^m,nV 


T^  n+,  + SIm-DAz  T„.,  „ + 
L Az  2Ar 

+ OfSImlAz  + p - h^(m)Az  + T^pjpcA|m)A 
k 2k  ( 2Atk 


pcA(m)Az 


hS(m)AzT. 


pcA(m)Az  — A(m) 


— S{m*1)Az  — hS(m)Az 


Cm  7 m - MR.  1 < n < MHL 


Thii  caw  covan  the  bottom  aurfeoa  of  tha 

avadge.  Tha  haat  balanca  li 

< 


Atk 


i 


'*"m,n-1  '*■  "^m,n+1  ■*■  Sfm-DAz  Tp,.i  „ 

Az 


pcA(m|Az  L ^ ^ Ar 

+ hS(m)AzTa  + ^m)A2  + p - h^SMAz  + pcA(m)Az 


k k k 

- 2A(m)  — S(m-1)Az  — hS(m)Az\*| 


V 


Az 


Ar 


7J 


Atk 


(11) 


Caw8  m>MR«n»MHL 


A(m) 


This  final  case  covers  the  right  lower  comer 
node,  the  only  node  with  two  outside 
surfeces.  The  heat  balance  is 


^"^m,n^new 


2Atk  Aim)  T,„  p.i  + S(m-1)Az  T^.^  „ + h |(A(m)  + S(m)Az)l  T, 
/)CA(m)A2  Az  5Ar  If  L o J 


+ q^^lm)  + S(m)  A^  + p - h,  |A(m)  + S(m)  AzJ-*-  T^ 

_h  |A(m)  + S(m) 


pcA(m)Az 


- A(m)  - S(m-1)Az  — 
Az  2Ar 


2Atk 


(12) 
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The  choice  of  numerical  velues  for  Ar,  Az  end  At  it  made  according  to  the 
convergence  criteria  of  Holmen  <see  reference  3).  The  coefficient  of  T^  ^ in  equation  {3) 
should  be  equal  to  or  greater  than  zero.  This  gives 


pc  2 2 

>0  (13) 

kAt  Ar*  Az* 

Any  combination  of  At,  Az  end  Ar  which  satisfies  (13)  will  insure  convergence  of  the 
numerical  solution. 

III.  MICROWAVE  HEAT  ABSORPTION 

Due  to  the  exponentiel  decay  of  microwave  energy  as  it  penetrates  the  roast,  each 
node  absorbs  a different  amount  of  energy.  In  tha  ensuing  discussion  microwave  energy 
will  be  discussed  in  terms  of  individuel  "beams"  of  energy.  In  any  microwave  oven  multiple 
reflections  of  energy  from  oven  wells  will  cause  incident  microwave  beams  to  strike  the 
roast  surface  with  angles  from  the  surface  normal  varying  between  0 and  90**.  Microwave 
beams  that  are  incident  et  any  angle  greeter  than  0°  ere  strongly  refracted  since  the  eir/meat 
interface  has  an  index  of  refraction  of  about  7 as  seen  in  refererwe  1.  Thus  a beam 
with  en  incident  angle  of  70°  would  actually  penetrate  the  meat  et  about  8°  from  the 
surface  normal.  Consider  a roast  cross  section  with  only  5 radiel  nodes  as  shown  in 
Figure  2 • 

A circle  is  drawn  halfway  between  each  node  defining  e "ring"  that  is  associated 
with  eech  node.  To  define  the  microwave  power  of  incident  beam  Pq  one  must  consider 
first  the  total  power,  Pt#  ebsorbed  by  the  roast.  This  power  is  assumed  to  be  equally 
distributed  over  the  roast  surface.  If  the  cross  section  shown  in  Figure  2 is  essumed 
to  be  a disc  with  thickness  Az,  the  fraction  of  totel  incident  power  received  by  this 
disc,  denoted  by  RT,  is 


RT  = r ^■*-1- - — .1  r —1  = 

L2irrl  + 2irr*  J L • J l+r 
Then  the  microwave  power  incident  upon  the  disc  would  be 
Pq  ’ <Pt)(RT)  cal/s 


(14) 


(15) 


As  seen  in  Figure  2,  en  incident  beam  with  power  P^  is  refracted  end  travels  obliquely 
through  rings  5,  4,  3 end  2;  it  does  not  trevel  through  ring  1 end  for  this  specific  beem 
ring  1 does  not  receive  eny  energy.  Microweve  powar  ebsorbed  is  calculated  according  to 
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where  P * remaining  power  at  depth  d (cm) 
pQ  • power  transmitted  through  surface 


(16) 


a ~ dielectric  attenuation  coefficient 


The  calculation  of  a is  discussed  in  reference  1. 


REFRACTED 


mm 


SUAFACf 


NDRMAi 


INCIDENT 


MICROWAVE 
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In  the  computer  progrem,  prior  to  eny  power  ibiorption  caiculetions,  the  velue  of 
co^  for  each  rh^  is  calculeted  as  riiown  in  Figure  3.  To  calculate  power  ebscMted,  one 
must  start  at  the  surface  (outermost  ring)  using  equation  (16).  Then 

P,  - P<,e"2aDs 
end 

« p^g— 2aD4  ■ « p^ef~2aD4~2aDjI 

accordingly 

Pj  . p^e-2laD3+aD4+0fDjJ  (17) 

Equation  17  cannot  be  further  simplified  as  a is  e function  of  temperature,  and 
each  ring  (or  each  node)  is  et  e different  temperature.  The  microwave  power  absorbed 
in  the  first  pessthrouflh  ring  3,  APj,  is 

AP,  - p4  - p,  * Pq  (6^2(004*^05)  _6-2(a03+a04+o0j)]  (,qj 

In  a simitar  manner,  the  power  absorbed  in  the  second  pass  through  ring  3,  AP'3 
is  calculated  as 

m p'j  ..  P'^  ■ p^  jg—2(aO5+aO4*HJtO3+20£OlS)_^— 2(005+004+003 +2aOlS)j 

(19) 

Then  the  total  energy  absorbed  in  ring  3 is  the  sum  of  energy  absorbed  in  the  two  passes, 
and  is  written  as 

P,j  - APj  + AP',  (20) 

The  power  ebeorbed  by  each  ring  can  be  calculeted  in  e similer  manner.  Referring 
to  Figure  1,  it  is  seen  that  only  e fraction  of  the  power  absorbed  in  each  ring  is  absorbed 
by  the  wedge^haped  nodal  area.  This  fraction  is  the  wedge  outer  arc  length  (assumed 
to  be  Ar)  divided  by  the  roast  outer  circumference  (here  SvAr,  see  Figure  2),  or  1/8ir. 
Then  the  radial  component  of  power  absorbed  by  node  3 would  be 

^node  3 “ ca'/*  (21) 

In  actual  operation  the  POWER  subroutine  calculetes  power  absorbed  in  e similer 
manner  but  somewhat  modified  to  make  computations  in  the  computer  more  efficient. 
To  gat  good  computed/experimental  deta  agreement  several  beam  incidence  angles  (^|) 
were  assumed  with  e certain  percentage  of  total  power  associated  with  each.  I n reference  1 
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a typical  sat  of  with  associated  percentages,  called  the  m/w  distribution,  was  40%  power 
at  01  “0®,  30%  power  at  0i  » 10®,  and  30%  power  at  03  » 20®.  For  each  0|,  an 
array  COBli,  J)  is  computed  (where  J » 1,MR)  which  stores  the  value  of  co0  for  each 
ring  (J)  for  each  incidence  angle  (I).  > 

Since  the  heat  transfer  portion  of  the  program, subroutine  INTER,  begins  with  node 
(1,1)  and  increments  radially  out  to  node  (MR,1),  then  node  (1,2)  to  (MR,2),  etc.,  the 
POWER  subroutine  first  calculates  the  power  remaining  at  the  most  interior  node  tta 
beam  passes  through  and  works  its  way  to  the  surface  by  doing  things  in  reverse. 

Microwave  power  entering  the  roast  ends  is  a much  simpler  case  to  handle.  A 
procedure  utilizing  a direct  application  of  equetion  (16)  is  used  for  energy  propagating 
in  the  z-direction.  Again  microwave  beams  are  incident  in  a random  manner,  with  angles 
from  the  surface  normal  varying  between  0®  ar>d  90®.  Due  to  the  large  refraction  at 
the  roast  surhK;e  end  the  plane  surface  involved,  it  is  assumed  that  the  microwave  erwrgy 
is  in  the  form  of  a plane  wave.  Then  the  simplifying  assumption  that  all  microwave 
power  is  incident  perpendiculetly  to  the  surface  can  be  made.  Energy  entering  both  ends 
is  taken  into  account.  In  a procedure  similar  to  that  already  discussed,  the  program  first 
steps  through  equation  (16)  from  the  surface  to  tha  center,  and  works  through  the  iteration 
process  in  reverse. 

IV.  SURFACE  PHENOMENA 

A.  Convective  Heat  Transfer 

The  temperature  difference  between  the  roast  surface  and  surrounding  air  causes 
tuitural  convection  currents  to  transfer  some  heat  from  (or  to)  the  roast.  The  complex 

process  of  convective  heat  transfer  is  summed  up  in  one  term,  h,  the  convectiva  heat  | 

transfer  coefficient.  The  measurement  of  h is  difficult  as  both  h and  h^  (evaporative  i 

cooling  coefficient,  discussed  in  the  next  section)  are  surface  cooling  phenomena  present  I 

during  mkrowave  cooking.  Early  research  on  meat  roasting  in  a heated  oven  (no 

microwaves)  with  the  roast  suspended  from  a scale  provided  evaporative  cooling  data, 

typically  as  shown  in  Figure  4.  Using  this  experimental  data,  a trial*and-error  process 

was  used  to  determine  h,  by  matching  the  computed  and  experimental  time-temperature 

profiles.  For  the  case  shown  in  Figure  4,  the  corresponding  temperatu’e  profila  is  shown 

in  Figure  5.  The  trial-and-error  process  gave  the  value  of  h to  be  0.(XX)11  C8i/g®Ccm^. 

This  value  of  h was  confirmed  by  several  more  tests,  and  was  used  throughout  all  research. 

B.  Evaporative  CooNng 

The  evaporative  cooling  coefficient,  h^  (called  HEVAP  in  computer  program)  has 
been  a difficult  paran>eter  to  work  with.  Without  the  benefit  of  surface  evaporation  and 


EVAPORATIVE  LOSS,  GRAMS  H^O 


COOKING  TIME,  MINUTES 


Flgiir*  4.  Evaporativ*  Lom  in  121"C  Ovan,  Baaf  RoaM  3380  9,  Surfaoa  Arm 
1186  cm^ 


L ^ 


ROAST  TEMPERATURE.  ®C 


COOKING  TIME.  MINUTES 

Figui*  5.  RoMt  Tamptritura  Histofy  in  121^*0  Ovm,  Batf  RoKt  3380  g 
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its  cooling  effect,  the  meet  neer  the  surfece  would  overcook  resulting  in  toughening, 
hardening,  and  unpalatability.  Acquisition  of  evaporative  cooling  data  during  microwave 
cooking  was  not  possible  as  the  roast  had  to  be  rotated  during  microwave  heating  and 
could  not  be  suspended  from  a scale.  U^ng  h ■ 0.00011  cal/g^Ccm^  as  previously 
discussed,  a trial-and-error  procedure  was  used  to  determine  the  value  of  h^  again  by 
matching  the  experimental  and  computed  temperature  profiles.  For  several  roasting 
conditions,  h^  is  shown  as  a function  of  surface  temperature  in  Figure  6.  The  data  of 
Figure  6 was  graphed  to  show  both  slope  and  constant  value  of  HEVAP  v.s.  wattage, 
providing  a linear  relationship  between  these  variables.  This  information  was  directly 
programmed  In  the  EVAP  subroutine  so  that  proper  velues  of  the  slope  and  the  constant 
value  of  HEVAP  are  automatically  assigned  given  the  frequency  and  wattage  of  microwave 
power  used.  For  a heated  oven  without  use  of  microwave  power,  or  combinations  of 
microwaw  heating  and  radiant  heating,  values  of  variables  HP  end  HMX  ara  read  in  as 
input  data,  and  HEVAP  « HP  • T0iM,N)  for  roast  surface  temperatures  less  than  40°C, 
and  HEV^  • HMX  for  surface  temperatures  greater  than  40‘’C. 

C.  Radiant  Heat  Transfer 

A detailed  analysis  of  radiant  energy  transfer  was  performed.  Four  surfeces  were 
involved:  (1)  the  roast,  <2)  the  oven  ends,  (3)  the  heating  elements,  and  <4)  the  curved 
oven  well.  The  size  and  shape  of  each  surfece  was  measured,  and  shape  factors  between 
all  surfaces  were  calculated.  Then  a set  of  4 equations  in  4 unknowro,  the  radiosities 
J|,  were  written  for  the  measured  surface  temperatures  for  each  oven  setting.  Using  the 
calculated  radiosities,  radiant  heat  transfer  to  the  roast  was  calculated  by  multipiving  the 
difference  in  radiosities  by  the  shape  factor.  Appendix  A shows  the  details,  and  the 
slopes  and  intercepts  of  the  curves  of  Figure  A3  are  included  In  element  MAIN  2. 

D.  Index  of  Refraction 

Inddent  microwave  radiation  nrperiences  a large  refraction  as  it  enters  the  roast 
surface.  Assuming  meet  has  a negligible  magnetic  loss  and  a permeability  equivalent  to 
that  of  a vacuum.  Von  Hippel*  givas  the  index  of  refraction,  n^,  as 


*Von  Hippel,  A.  R.,  Dialectries  and  Wavas,  MIT  Press,  Cambridge,  MA.  1954. 
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Figur*  0. 


Evqionttiva  Cooling  Cosfficiont  Dotarminod  by  Computad/Expari- 
mantal  Data  Matching,  with  h ■ 0.00011  c8l/cm^t°C 


23 


where  > dielectric  constant  ' (22) 

f"r  dielectric  loss  fector 

For  fresh  raw  beef  the  velue  of  n,.  over  the  range  of  0— 100°C  ranges  from  6.5  to 
7.7  for  frequencies  from  900  to  2800  MHz.  For  frozen  beef,  the  value  of  n^  is  2.1 
to  2.2  for  ^ above  frequency  range.  For  example,  if  a microwave  "beam"  was  incident 
at  46“  from  the  roast  surface  normel,  it  would  penetrate  the  meat  at  approximately  19“ 
for  frozen  meat,  and  approximately  6“  for  thawed  meet.  Thus  the  microwave  distribution 
end  an>ount  of  focusing  of  microwaves  is  much  less  for  a frozen  roast  than  e thawed 
one.  It  is  important  to  note  that  the  index  of  refraction  depends  on  the  surface 
temperature,  and  is  an  important  consideration  when  using  microwave  energy  to  thaw 
roasts. 

V.  SUPPORTING  DATA 


A.  Dielectric  Properties 


The  dielectric  properties  of  meat  differ  for  the  two  microwave  frequencies  used,  915 
and  2450  MHz.  From  the  values  of  cV  end  e'V  reported  by  Goldblith  end  Weng*  for 
915  MHz,  the  velue  of  a,  the  attenuetion  coefficient,  was  calculated  according  to 


(23) 


where  v » microwave  frequency,  cycles/s 


2irv  ^ T 

T 


fl  + e'V/eV 


-j 


Cy  ~ speed  of  light  in  vacuum,  cm/s 

The  value  of  a wes  plotted  as  a function  of  temperature  (see  reference  1),  divided  up 
into  10“C  segments,  end  a median  value  of  a assumed  constant  over  each  10“C  interval. 
Thus  the  variation  of  a with  tempereture  is  approximated  by  a stepped  curve;  this 
information  is  staed  in  an  array  in  computer  function  ALPHA.  When  the  value  of  a 


! 


i 


* Goldblith,  S.  A.,  aiHl  War>g,  D.I.C.,  "Dielectric  Properties  of  Foods,"  US  Army  Netick 
R4D  C^mand,  Technical  Report  TR  76— 27-FEL,  1975. 
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for  a particuler  node  is  required  in  computer  computations,  function  ALPHA  is  called. 
The  node  temperature  is  divided  by  10,  this  result  truncated  giving  en  integer,  end  the 
value  of  01  associated  with  that  integer  is  used. 

For  2450  MHz,  dielectric  data  from  Sengtsson  and  Risman‘  and  Ohisson  and 
Bengtsson^  were  used.  Although  obtained  for  2800  MHz,  the  data  were  used  for  2450 
MHz  (see  reference  1 for  the  explanation).  Equation  (23)  is  used  to  obtain  values  of 
a,  and  the  same  step-epproximation  previously  described  is  used  tn  digitize  the  data.  Here 
function  ALPHF  is  used  to  store  the  data  and  obtain  the  curr‘*<  t value  of  a (according 
to  the  node  temperature)  when  called. 

B.  Thermal  Conductivity 

The  variation  of  the  thermal  conductivity  of  meat,  k,  with  temperature  and  percent 
moisture  was  accounted  for  in  this  research.  Data  from  Hill*  was  linearized  with  respect 
to  percent  moisture,  and  nonlinear  equations  were  derived  for  k as  a function  of  both 
temperature  and  percent  moisture.  Figure  7 shows  the  raw  data  plotted  and  Figure  8 
shows  the  variation  of  the  slope  and  intercept  of  the  Figure  7 curves  plotted  as  a function 
of  temperature.  Thus  for  a given  temperature,  from  the  equations  of  the  Figure  8 a 
value  of  slope  and  intercept  is  found  which  in  essence  defines  a straight  line  on  Figure  7; 
then  given  the  percent  moisture,  a unique  value  of  k is  identified  in  Figure  7.  Subroutine 
COND  contains  the  equations  illustrated  in  Figures  7 and  8. 


*Bengtsson,  N.  E.,  and  Risman,  P.  O.,  "Dielectric  Properties  of  Foods  at  3 GHZ  as 
Determined  by  a Cavity  Perturbation  Technique  II  Measurements  on  Food  Meteriels,"  J. 
Microwave  Power,  6 (2),  1971. 

’Ohlson,  T.,  and  Bengtsson,  N.  E.,  "Dielectric  Food  Data  for  Microwave  Sterilization 
Processing,"  J.  Microwave  Power,  10  (1),  1975. 

*Hill,  J.  E.,  Leitmen,  J.  D.,  and  Surxlerland,  J.  E.,  "Thermal  Conductivity  of  Various 
Meets,"  Food  Tech.,  21  (8),  1967. 
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Figure  8.  Fitted  Curvet  for  Slope  and  Intercepts  from  Figure  7 vs  Roast 
Temparature 


27 


r 


C.  Enttialpy  and  Spaeific  Haat 

Enthalpy  is  a term  used  to  describe  the  change  in  energy  required  to  heat  meat 
from  a frozen  state  to  a thawed  state  (meat  is  over  70%  water)  ^ich  occurs  ovar  a 
tamperatura  range  of  about  -10^  to  0”C.  Dickerson’  determined  values  of  enthalpy  (units 
cal/g)  for  meat  over  the  temperature  range  -4(f  to  The  instantaneous  slope  of 

this  curve  is  the  valua  of  specific  heat,  c,  (units  cal/g^C)  at  the  corresponding  temperature. 
The  enthalpy  curve,  shown  in  Figure  9,  was  linaarized  over  9 temperature  intervals,  each 
interval  with  a different  slope,  or  specific  heat.  This  temperature  and  slope  data  is 
incorporated  into  subroutine  SPHT.  This  subroutine,  when  called,  identifies  the  proper 
interval  for  the  given  node  temperature  and  chooses  the  proper  value  for  the  specific 
heat. 

Vt.  COMPUTER  PROGRAM  SUMMARY 
A.  Elements 

The  program  file  is  named  MEAT* ROAST,  and  consists  of  a main  program  and  8 
sub-programs.  The  elements  are  described  as  follows: 

Element  Subroutine/Function 


Name 

Name 

Description/Function 

MAIN2 

Main  program  which  calls  all  subprograms, 
carries  out  iterative  process,  reads  input  data, 
prints  output. 

MAIN1 

Main  program  which  is  used  to  obtain 
^dimensional  plots  - same  MAIN2  except  in 
way  output  data  handled. 

INTER 

INTR 

Calculates  intarior  node  temperatures  using 
finite  difference  approximation. 

BNDRY 

BNDR 

Calculates  boundary  node  temperatures  taking 
into  account  radiant  and  microwave  energy. 

’Dickarson,  R.  W.,  Jr.,"Tharmal  Properties  of  Foods,"  in  The  Freezing  Preservetion  of 
Foods,  4th  edition,  Avi  Publishing  Co.,  Westport,  CT,  1968. 
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B.  input  Data  Format 


There  ere  10  lines  of  input  data  necessary  to  completely  define  roasting  conditions. 
These  10  lines  ere  contained  in  a data  file  called  DATA,  and  are  in  the  free  format; 

i.e.,  reai  numbers  separated  by  a comma.  DATA  variabies  must  be  in  the  following  order: 

1.  P2450,  P915,  TIN,  TAM,  ENDPCT,  NSHLD 

2.  D1,  D2,  XL,  AZ,  AT.  TMAX,  KDELO 

3.  RHO,  H,  PH20 

4.  TAM1,  TAM2 

6.  T2450A,  P2450A,  T24508.  P2450B,  T9 15,  PN915,  TOVEN.  TMPOVN 

6.  HP,  HMX 

7.  NBT 

8.  TT(I),  PCTT(I) 


9.  NBF 


j 

i- 

t 

I 


10.  TF{I).  PCTFd) 


Definition  of  the  above  variables  is  as  follows: 


Variable 

P2450 

P915 

TIN 

TAM 

ENDPCT 

NSHLD 

01 

D2 

XL 

AZ 

AT 

TMAX 

KDELO 

RHO 

H 

PH20 
TAM1  ' 
TAM2^ 


Definition 

Oven  power  at  2450  MHz,  watts 
Oven  power  at  915  MHz,  watts 
Initial  temperature  of  roast, 

Oven  ambient  temperature,  "C 

Fraction  of  total  microwave  energy  entering  roast  through  shield 

Number  of  node  where  shield  overlap  begins 

Diameter  of  raw  roast,  cm 

Diameter  of  cooked  roast,  cm 

Roast  length,  cm 

Axial  node  spacing,  cm 

Time  increment,  s 

Maximum  cooking  time,  s 

Printout  interval  in  number  of  time  increments 

Meat  density,  g/cm^ 

Corwective  heat  transfer  coefficient,  cal/s^’C  cm^ 

Percent  moisture  in  meat 

Used  to  compute  oven  temperature  rise  (about  20—50°  C) 

When  cooking  with  microwaves,  TAM  * TAM2  + TAMl  * TMIN 


] 

j 
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Varilbto 


Dafinition 


T24S0A  Tima  at  which  P2450  h to  change,  min 

P2450A  New  valua  of  P2450,  watts,  at  T2460A  min 

T24S0B  Tima  at  which  P2460A  it  to  change,  min 


P2450B 


New  valua  of  P2460B,  watts,  at  T2450B  min 


Tima  at  which  P915  is  to  change,  min 


PN915 


New  valua  of  P915,  watts,  at  T915  min 


TOVEN 


Tima  at  which  TAM  is  to  change,  min 


TMPOVN  New  value  of  TAM,  ‘‘C,  at  TOVEN  min 

HP  Up  to  40®C  roast  surface  temperatura,  HEVAP  - HP  * TO(M,N) 

HMX  Greater  than  40*’C  roast  surface  temperatura,  HEVAP  » HMX 


PCTTII) 


Number  of  incident  angles  of  incoming  radiatioifj 
Angle  of  incoming  radiation  to  surface  normal,  rad.  \thawed 


Percent  of  total  radiation  at  angle  TT(I) 


PCTF(I) 


Same  as  NBT,  TT(I)  and  PCTT(I)  except  for  frozen  roast 


Several  features  of  the  program  not  previously  mentioned  can  be  convaniantly 
explained  now. 

(1)  Partial  or  full  shielding  of  tha  roast  ends  is  possible  using  the  ENDPCT  input. 
Shield  ovarlap  can  also  be  specified  using  NSHLD. 


(2)  Radial  shrinkage  of  tha  roast  during  cooking  is  accounted  for  by  inputs  D1 
and  D2.  Twenty-four  equal  diametral  divisions  are  assumed  and  linear  shrinkage  from 
diameter  D1  at  start  to  diameter  D2  at  finish  is  accounted  for  by  modifying  Ar  each 
iteration. 
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(3)  Chengos  in  microwave  power  or  oven  temperature  can  be  made  et  e specified 
time  using  input  line  5. 

(4)  If  the  roast  is  initially  frozen  (i.e.^.TIN  < 0*’C)«  proper  choice  of  property 
variables  and  microwave  distribution  is  made  automatically  by  IF  statements  within  the 
program. 

(5)  An  optionel  30  minutes  cooling  (thermal  equilibrium)  period  has  been  provided 
for  in  MAIN2.  During  this  time  the  roast  is  outside  the  oven«  the  temperature  profile 
tends  to  even  out.  An  outright  guess  for  values  of  h end  h^  was  used;  for  relieble  results, 
values  for  h end  he  ^ ^ determined  experimentally. 

C.  ThrM-Oimensional  Plotting 

The  CALCOMP  plotter  is  used  to  plot  a threeKlimensional  surface  representing 
temperatures  over  e center  exiel  cross-section  of  the  entire  roast.  After  performing  all 
the  operations  contained  in  MAIN2,  elentent  MAIN1  manipulates  the  temperature  data 
for  the  1/4  cross-section,  using  symmetry  considerations  to  create  a temporary  file 
(called  15.)  with  temperatures  representative  of  e full  axial  cross-section.  TheTHREE  D/ll 
Calcomp  software  system  is  used;  this  system  requires,  for  our  case,  e data  input  of  1 1 lines 
which  specify  ell  characteristics  of  the  final  3-0  graph.  The  MAIN1  program  is  set  up 
to  store  temperature  profile  at  preset  intervels,  typically  5 minutes,  in  the  file  15.  When 
the  roast  is  finished  cooking  (mathematically),  file  15.  is  quite  large,  containing  typically 
10  to  20  sets  of  temperature  profiles  to  be  plotted.  In  the  THREE  D/II  data  input, 
line  A tells  the  plotter  how  many  plots  to  expect  (along  with  other  information)  end 
Lines  B through  M contain  graph  characteristics.  Thus  two  elements  need  to  be  used 
to  edd  the  THREE-D/II  input  data.  Element  ROAST.B-M  contains  10  lines  of  input 
for  lines  B through  M;  element  ROAST.ADDS  contains  line  A input,  followed  by 
statement  ADD  ROAST.B-M"  10  to  20  times,  one  time  for  each  plot  desired.  Thus 
the  THREE-D/II  software  will  produce  e three-dimensionel  temperature  plot  every  preset 
interval  of  cooking  time.  An  example  plot  is  shown  in  Figure  10.  An  animated  film 
showing  sequential  temperature  profile  plots,  from  start  to  finish  of  cooking,  has  been 
made  for  several  heating  methods.  This  film  shows  very  dearly  the  difference  between 
heating  methods  throughout  the  entire  roast. 

VII.  SUMMARY 

A versatile  heet  transfer  program  for  computer  simulation  of  redient/convective 
Iweting  of  cylirxJrical  food  substances  has  been  presented.  Provisions  for  microwave  end 
shielding,  radial  shrinkage,  changes  in  heating  method  during  cooking,  thawing,  and  e 
temperature  equilibration  period  have  been  included.  Detailed  heat  transfer  equations 
and  computationel  techniques  have  been  shown. 
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Figur*  10.  ThrM-DinMMional  Tempareture  Plot  for  Microwsva  Hotting  «t  300  Watts 


at  2450  MHz,  60  Mimitos  Hooting  Timo 
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APPENDIX  A 

RADIATION  HEAT  TRANSFER 

In  order  to  carry  out  a radiation  heat  transfer  analysis,  the  shape  factors  for  tha 
various  heat  transfer  surfaces  need  to  be  computed.  These  are  obtained  from  a 
consideration  of  the  oven  and  roast  geometry,  as  shown  in  Figure  A-1: 


A1:  ROAST 

A2:  OVEN  ENOS 

A3;  HEATING  ELEMENTS 

A4:  OVEN  CURVED  SURFACE  (HEATED) 


0.117  m2 
0.280  m2  EACH 
0.155  m2  each 
1.452  m2 


Figure  A1.  Oven  cavity,  59.7  cm  dia  by  94  cm  long,  showing  heat  transfer  surface 

8f^O0S« 
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AI.Kreith,  F.,  Prindfyhs  of  Heat  Transfer,  international  Textbook  Co.,  Scranton,  PA,  1967. 
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In  the  electrical  analogy,  this  equation  is  interpreted  as  a rate  of  current  flow  between 
nodes  Ej^i  and  Ji  connected  by  a resistance  Pi/Aiej.  Connecting  nodes  J]  through 
J4  by  resistances  due  to  shape  factors  gives  the  electrical  analogy  shown  in  Figure  A2: 


1 

^2^2  Jz  *2^23  J3  ^^3 


Figure  A2.  Electrical  Anirtogy 

Here  Ej^  is  the  black  body  emissive  power  of  surface  1 and  is  equal  to  oTi  * where 
T|  is  the  at^lute  temperature  of  surface  1 and  o i»the  Stefen-Boltzmann  constant.  The 
radiosity,  Ji,  is  defined  as  the  net  rate  at  which  radiation  leavessurface  1 per  unit  area. 


To  reduce  the  electrical  analogy  network  to  a set  of  equations  to  solve,  an  equation 
balancing  the  currents  at  each  radiosity  node  is  written.  For  node  Ji  this  is 


I 


This  reduces  to 

1.11J,  - .026Jj  - .OI4J3  - .O68J4  - E. 


(3A) 


Similar  equations  written  for  die  other  nodes  Jj,  Js,  and  J4  can  be  combined  with 
equation  (3A)  into  matrix  form,  as 


1.11 

-.026 

-.014 

-.068 

-.046 

1.9S 

-.14 

-.77 

-.005 

-.026 

1.10 

-.067 

-.046 

-.29 

-.133 

1.47 

Jj 

J4 


■b4 


(4A) 


Then  for  known  values  of  through  Ej,  , which  are  calculated  from  the  surface 
temperatures,  solution  of  the  four  simultaneous  equations  yields  values  of  the  radiosities. 
Then  the  radiant  heat  transfer  received  by  the  roast  is  calculated  by 


— roast  “ (Jj  ~ J| ) F|  .1  + (Jj  ~ Ji  IF|  *$  + (J4  ~ .4 

A 


(5A) 


Surface  temperatures  T3,  T^,  and  T4  for  various  oven  temperature  settings  ware 
determined  by  experiment  Table  A1  depicts  typical  values  obtained: 


OVEN 

SETTING 

Tj 

T, 

T4 

'C 

K 

K 

K 

K 

121 

394 

389 

458 

333 

Table  A1. 

149 

422 

417 

497 

355 

Surface  Temperatures 

177 

450 

444 

550 

378 

for  Various  Oven 

204 

478 

478 

611 

403 

Settings 

41 


i 


Three  repreientativa  roast  temperatures  are  chosen,  15.6**C,  54.4'*C,  and  93.3° C.  For 
these  roast  temperatures,  a value  of  radiosity  for  each  ovan  setting  is  found  by  solving 
the  sat  of  equations  {AA).  Table  A2  summarizes  the  results. 


Oven  Temperature 
Setting  (°C) 


Radiosity 


cal 


Scm^ 


Node 


For  Roast  Temperature  (°C) 
15.6  54.4  93,3 


121 

.0117 

.0174 

.0246 

^2 

.0254 

.0257 

.0260 

Js 

.0570 

.0570 

.0571 

J4 

.0316 

.0319 

.0322 

149 

Jl 

.0129 

.0185 

.0250 

J2 

.0336 

.0338 

.0342 

J3 

.0787 

.0788 

.0789 

J4 

.0420 

.0422 

.0425 

177 

Jl 

.0145 

.0200 

.0281 

.0450 

.0452 

.0455 

Jj 

.1173 

.1173 

.1174 

J4 

.0559 

.0561 

.0664 

204 

Jl 

.0169 

.0225 

.0305 

J2 

.0619 

.0621 

.0624 

J3 

.1781 

.1781 

.1782 

J4 

.0769 

.0771 

.0775 

Table  A2.  Radiosity  Values  for  Various  Oven  Temperature  Settings 
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A,-A4 


NOMENCLATURE 


Area  associated  with  node  cm^ 

Meet  transfer  surfece  areas' in  oven, 

Specific  heat,  cal/g^C 

Speed  of  light  in  e vacuum,  m/s 

Depth  measured  from  meat  surface,  cm 

Distance  microwave  beam  travels  through  ring  m,  cm 

Black  body  emmissive  power  of  surfaces  1 through  4 
Shepe  factor  for  surface  1 to  surface  4 
Convective  heat  transfer  coefficient,  cal/s‘C  cm^ 

Evaporative  cooling  coefficient,  cal/s  cm^ 

Radiosities  of  surfaces  1 through  4 
Thermel  conductivity,  cal/cm  s^C 
Redial  node  component 
Maximum  numerical  value  of  redial  node 
Maximum  numerical  value  of  axial  node 
Axial  node  component 
Index  of  refraction 

Microwave  power  absorbed  in  volume  associated  with  node  (m.n),  cal/s 

Microwave  power  remaining  at  depth  d in  roast,  cal/s 

Microwave  power  transmitted  through  roast  surface,  cal/s 

Microwave  power  remaining  in  beam  efter  first  pass  through  ring  m,  cal/s 

Microwave  power  remaining  in  beam  after  secorKt  pass  through  ring  m, 
cat/s 
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1 


1 


^mT 

Pt 


« 

<1 


f 


r 


RT 


S(m) 

S(m-1) 

t 

T 

T,~T4 

T, 

V 


z 

a 

? 

A 


♦ 

0 

fi 


P1-P4 


Total  microwave  power  ioit  in  ring  m (any  integer),  cal/t 
Microwave  power  abeorbed  by  entire  roait,  cai/s 
Rate  of  heat  addition,  cai/s  cm’ 

9 

Rate  of  redient  heat  addition,  cai/s  cm’ 

Radial  axis 

Fraction  of  total  power  received  by  disc  with  thickness  Az 
Outer  arc  length  of  area  surrounding  node  (m,n),  cm 
inner  ere  length  of  area  urroundi^  node  (m,n),  cm 
Time,  s 
Temperature, 

Absolute  temperatures  of  surfaces  1 through  4,  K 
Ambient  temperature  urrourxiing  roast,  °C 
Temperature  of  node  (m,n), 

Volume  asiocieted  with  noda  (m,n),  cm’ 

Axial  axis 

Attenuation  coefficient,  cm*^ 

Angle  the  microwave  beam  makes  with  radial  line  in  given  ring 

Increment;  used  as  e prefix,  e.g.,  Ar 

Refracted  beam  angle,  degrees  from  normal 

Angular  axis 

Meat  density,  g^cm’ 

Reflectivity  of  surfaces  1 through  4 
Dielectric  constant 
Dielectric  loss  factor 
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tt’~€4  Eminivity  of  lurfiOM  1 through  4 

V Microwav*  ffaquancy,  cyclei/s 

o Sttphan  — Boltzmann  constant 
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